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At very low temperatures, the tunneling theory for amorphous solids predicts a thermal conduc- 
tivity K <x T^, with p = 2. We have studied the effect of the Nuclear Quadrupole moment on the 
thermal conductivity of glasses at very low temperatures. We developed a theory that couples the 
tunneling motion to the nuclear quadrupoles moment in order to evaluate the thermal conductivity. 
Our result suggests a cross over between two different regimes at the temperature close to the nu- 
clear quadrupoles energy. Below this temperature we have shown that the thermal conductivity is 
larger than the standard tunneling result and therefore we have p < 2. However, for temperatures 
higher than the nuclear quadrupoles energy, the result of standard tunneling model has been found. 

PACS numbers: 61.43.Fs, 65.60.-|-a 



I. INTRODUCTION 

Amorphous or glassy materials differ significantly from 
crystals, especially in the low temperature range. Be- 
low IK, the specific heat Cv of dielectric glasses is much 
larger than in crystalline materials. Moreover, the ther- 
mal conductivity k is orders of magnitude lower than 
the corresponding values found in their crystalline coun- 
terparts. Cy depends approximately linearly and k al- 
most quadratically on temperaturei. The generally ac- 
cepted basis to describe the low temperature properties of 
glasses is the phenomenological tunneling model^i^. To 
explain these behaviors, it was considered that atoms, 
or groups of atoms, are tunneling between two equilib- 
rium positions, the two minima of a double well poten- 
tial (DWP). The model is known as the two level system 
(TLS). In the standard TLS model, these tunneling ex- 
citations are considered as independent, and some spe- 
cific assumptions are made regarding the parameters that 
characterize theroSi. 

The TLS can be excited from its ground state to the 
upper level therefore contributed to the heat capacity. 
TLSs can also scatter phonons and in this way decrease 
their mean free path and, correspondingly, the heat con- 
ductance. 

New interest in this problem was stimulated by several 
experimental results^'^'^'^. Until these experiments it was 
the general believe that the dielectric properties of insu- 
lating non-magnetic glasses are independent of external 
magnetic field. It is very surprising that strong magnetic 
field effects were discovered in polarization echo experi- 
ments at radio frequency and in low frequency dielectric 
susceptibility measurement at very low temperatures^i^ii^. 
Several generalizations of the standard TLS model have 
been reported after the anomalous behavior of glasses 
in a magnetic field. According to these solutions, the 
models can be divided into " orbital" and " spin" 

models (nuclear quadrupole effect The "or- 
bital" models can provide an explanation for some of the 
magnetic field effects by considering the flux dependence 
of the tunneling splitting. Unfortunately, some assump- 



tions have been made which cannot be reconciled with 
the standard features of the tunneling model. 

A surprising outcome of these experiments is a novel 
isotope effect observed in different glasses^. The lat- 
ter effect shows the important influence of the nuclear 
quadrupole moments on the observed magnetic fleld de- 
pendence. Therefore it is very important to flnd the 
effect of nuclear quadrupole moments on the response 
function of glasses. For this purpose, in this paper we 
have studied the thermal properties of heat conductivity 
of cold glasses taking into account the quadrupole effects. 
In Section II using Wiirger's formalismi^, we introduce 
the nuclear spins in the frame of the two level system 
model. We will flnd the general form of the heat conduc- 
tivity of cold glasses which takes into account the nuclear 
quadrupole moment in Section III. And flnally in section 
IV, we end this paper by a summery and conclusion on 
our results. 



II. TLS COUPLED BY A NUCLEAR SPIN 

The standard TLS can be described as a particle or 
a small group of particles moving in an effective double- 
well potential. At very low temperatures only the ground 
states of each wells are relevant. Using a pseudo-spin 
representation the Hamiltonian of such a TLS read as 

Htls = ^Aocr^ + ^Aa,, (1) 

where A is the energy off-set at the bottom of the wells, 
and Ao is the tunnel matrix element. Diagonalization of 
this two state Hamiltonian gives the energies 

where E is the energy difference between the two wells. 
According to the randomness of the glassy structure, the 
energy difference between the two wells have a broad dis- 
tribution. The energy off-set and the tunneling matrix 
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element are widely distributed and are independent of 
each other with a uniform distribution of 
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where Pq is a constant. Using the notations u ■ 
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where {tp-) and \'4'+) are the ground and exited state of 
the system, respectively. For the moment there is no rig- 
orous theory for tunneling in glasses. It is assumed that 
atoms or groups of atoms participate in one TLS. As 
we mentioned before, in the case of the multi-component 
glasses, one or several of the tunneling atoms carry a nu- 
clear magnetic dipole and an electric quadrupole. When 
the system moves from one well to another, the atoms 
change their positions by a fraction of an Angstrom. 

We can describe the internal motion of the nuclei by a 
nuclear spin I of absolute value — h^I{I+l). For a nu- 
cleus with spin quantum number / > 1 the charge distri- 
bution p(r) is not isotropic. Beside the charge monopole, 
an electric quadrupole moment can be defined with re- 
spect to an axis e 



Q = J d^r [Sir- ef -r^] p{r). 
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Therefore each level of the pseudo spin projection will 
split to (21 + 1) nuclear spin projections with the quan- 
tization axis m = — /, ...,/. 

This can couple to an electric field gradient (EFG) at 
the nuclear position, expressed by the curvature of the 
crystal field potential. The potential describing this cou- 
pling is writteni^ 



Vq 



-eQ 
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■[Vnli+V22li + V33ll]. 
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The bases used here (61,62,63) are the principal axes 
of the tensor Vij which describes the electric field gra- 
dient, and 6 is the electron charge. According to the 
Laplace equation the potential obey V11 + V22 + V33 — 0. 
If we define the asymmetry parameter 77 — -^^^^j^, the 
quadrupole potential can be expressed as: 

Vq = e,[3I^ + ^{ll-ll)-l'] (6) 
'^'^^"s the quadrupole cou- 



4/(2/-l) 



where we denote by eg 
pling constant. 

Therefore we can write the quadrupole potential in 
terms of the reduced two-state coordinate: 



Ho = 



(7) 



where VJ^*-^^ is defined in Eq. ([5]) for the particles in 
right (left) welU^. We can go to basis \i^±[I ,'m±)) = 
|^±) ® \I:m±) which have defined as foUowing^^ 



H±\'ip±{I,m±)) = S±,m^|V'±(/, m±)) 
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where H± ~ 
fore E±^yn^ 
satisfy: 
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and there- 



±-j + em±; the corresponding eigen-states 
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and since i?+ and do not commute, their eigen-states 
are not generally orthogonal: 



(10) 



where these overlaps are dependent on the angle 9. (here 
6 is the angle between two axis of the Nuclear quadrupole 

in each wellsii: 6^,6^) 

III. HEAT CONDUCTIVITY 

The dominant effect of uniform strain field (describ- 
ing the interaction of the TLS with a phonon field) is 
on the energy of the tunneling state by changing the 
asymmetry energy. The changes in the barrier height 
can usually be ignoredi^. Any external perturbation is 
therefore diagonal in the local representation (|L), \R)) 
which when transformed into the diagonal representation 
( IV-"-)) has the form 

Hint = i^CTx + ^0-2)76 C0S(W<) = F-„t COs{ujt) (11) 

in the presence of a strain field ^ — ^ocos(ti;t), where 
^0 and uj are the amplitude and the frequency of the 
strain field respectively. The strain is given by e = ^q^qj 
and the parameter 7, defined as is equivalent to 

elastic dipole moment. Where is the phonon wave- 
vector with polarization a. Here the tensorial nature 
of 6 has been ignored and 76 is written as an average 
over orientations. Therefore we can easily sho'w^ that 
6 = (2^)^^"' where p is the bulk density and h is the 
Planck constant. 

Using Fermi Golden Rule, one can obtain the contri- 
bution of a phonon with wave vector and polarization 
a to the generalized TLS transition probability due to 
phonon emission and absorption, respectively: 

C>_„,_ - -J- I (V^+(/,mV) I H[^ I ^-(/,m_)) p 
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where ue 



e ^^i'-^i /Z is the Boltzmann weight, 



both wells are the same (xmi,m_ = ^m^.m- ^ ^m- 



constant and T is temperature. It must be noted here 
that the transition between the same TLS levels are zero: 



Kb is the Boltzmann ^m' ) , the famous result of the standard TLS model can 



(V._(/,m±)|ij;„JV+(/,O)=0 



r / = 



Therefore the phonon relaxation time can be found by 
summing — r^"° over all spin states: 
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pvl 
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where f q. is the sound velocity. Denoting 



(E) = nE_ 
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and after some calculations and averaging over TLS pa- 
rameters (using Eq. ^ , it can be easily shown that 



E 

m- .m'l 



I Xm^ ,m 
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Neglecting the phase difference between the nuclear mo- 
ments in the two wells and assumming that the EFG in 



be found: 



-1 _ PpTTjluJa 



(15) 



The thermal conductivity k{T) is evaluated on the as- 
sumption that heat is carried by non-dispersive sound 
waves, therefore one can write 

1 f°° 

k{T) = oE / l{'^a)Cv{l^a)g{u!a)VaduJa, (16) 
^ "'0 



where l{uJa) = Va/r^^ is the phonon mean free path 
of angular frequency to a, di^a) ~ 2^2° 3 is the phonon 
frequency distribution function, and Cv{<^a) is the heat 
capacity of phonon which is given by 



Cvil^a) 



1 ,hLj^ 
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and using the above equations the 



By defining x 
heat conductivity can be obtained, 

k{T) = S(T) X KTLs(r) (18) 

where ktls{T) = ^-^lp-^K%T^ is the standard 
TLS (STLS) heat conductivity^^, and the coefficient 
T.{T) is defined by 



4 f 



csch^{x)dx 
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As the exact behavior of the Heat Conductivity can 
not be found analytically, we are trying to solve Eq. [19] 
numerically. Assuming that I — 1 and — 1 mK as sug- 
gested by echo experiments, we observed the behavior of 
parameter S(r) in term of temperature. The results are 
peresented on Fig. ([T]) for different values of quadrupole 
angle {9) and by averaging over the 77 parameter. 

It can be seen that in high temperature regimes (e, <^ 
KbT), this ratio (E) goes to one. As it is predictable 
where the nuclear part effect can be neglected and heat 
conductivity behaves as the Standard TLS model. De- 
creasing temperature this ratio grows and will be satu- 
rated at very low temperatures. 

In agreement with expectation, at zero quadrupole an- 



gle the heat conductivity is the same as the result found 
from the standard TLS model (please see EqllJ] and the 
statements before that). At low temperature regime, in- 
creasing the quadrupole angle with small value cause the 
heat conductivity saturated value to be larger than what 
is found form the standard TLS model up to one percent; 

The same behavior can be found for / = 3/2 and 1 = 2. 
Also it can be shown that by changing eq, the growing 
regime shows a dependency on the quadrupole energy 
value; it means that by increasing the magnitude of e^, 
the growing regime will be shifted to higher tempera- 
tures. It shows that there is a cross over between two dif- 
ferent regimes in the temperature around the quadrupole 
energy value. 
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FIG. 1: Thermal variation of the ratio parameter E, for 
I = 1, eq — ImK and the different angle between the Nuclear 
Quadrupole moment in each TLS wells. 



IV. SUMMARY AND CONCLUSION 

In this paper we have studied the thermal properties of 
heat conductivity of cold glasses taking into account the 
quadrupole effects. To describe the interaction of a TLS 
with nuclear quadrupole, we have used a generalization 
of the standard TLS Hamiltonian with nuclear spin. The 
nuclear quadrupole of these systems leads to a splitting of 
the nuclear spin levels, which is different for the ground 
state and the excited state. The presence of this multi- 
level structure causes to increase the heat conductivity 
magnitude compared to that of simple two-level systems 
at lower temperature regimes {KbT ^ Cg). 

It is shown that the heat conductivity has a cross over 
between two different regimes at the temperature which 
is close to the nuclear quadrupole energy, KbT ~ e^. 
It can be observed that there are three different regimes 
versus temperature. 

The first regime deals with high temperature regime: 
KbT ^ Eq where the known standard tunneling model 
results have been found. At this regime the effect of 
quadrupole energy can be neglected in comparison with 
the TLS energy scale, therefore the nuclear spin splitting 
is not observable. 

The second regime demonstrates the temperature 
around nuclear quadrupole energy: KbT ~ e^- In this 
regime, by decreasing the temperature the heat conduc- 
tivity increases. This is the cross over regime between the 
standard TLS behavior and the low temperature regime 



where the nuclear quadrupole effects become important. 
In this area the nuclear quadrupole energy levels play 
an important role in the thermal behavior of the system. 
Decreasing temprature the nuclear quadrupole energy is 
comparable to thermal fluctuations. 

In general these sub-energy level are not the same 
in both wells of the TLS. Thus their eigen-states are 
not orthogonal and have the overlap with each other 
{\Xm'^,mA^ ^ Sm'_^,m_ and em+ 7^ em_)- This effect 
causes the mean free path of phonons to increase there- 
fore the thermal conductivity has larger value in compar- 
ison with the simple two level system at low temprature 
regime. This means that where KbT ~ e^, S = 1 -|- e{9) 
and heat conductivity exponent, p, is less than two in- 
stead of the p = 2 which has been found for standard 
TLS model. 

Finally for the third regime the heat conductivity will 
be saturated at KbT <c eq. 

To obtain a theorytical expression for this effect, one 



can write Xn 



and 



7m+ where e^^ 



em_ , and move to the special limitting 
case where <im+,m- < /37m+ , and 7^+ < e'^^ . By a httle 
manipulation it can be easily shown that 



e{e) = Cx {SPjf 



(20) 



where C = [AStt^ 
a numerical constant; 

-0e 



(/37) =E™^e---+/37™+ 



384C(3)] 



w 0.428 is 



-0< 



+ , and (/3^7^) = 

It shows clearly that 
by increasing the quadrupole angle the difference of sub- 
energy in both wells increases which causes the S(T) 
value to increase, in agreement with numerical results. 

In conclusion we believe that nuclear quadrupoles play 
an important role in the nature of glasses at low temper- 
atures. In this respect for solving the problem of cold 
glasses, it is useful to find the effect of nuclear spin on 
the other response functions. As far as we know there is 
no experimental result for the heat conductivity^! in the 
case KbT ^ eq. Therefore it might be a good suggestion 
for future experiments to approach lower temperatures 
or use the glasses with larger quadrupole energy. 
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